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Abstract 

We present the first lattice QCD calculation of the form factor for B — ► D*tP with three flavors 
of sea quarks. We use an improved staggered action for the light valence and sea quarks (the MILC 
configurations), and the Fermilab action for the heavy quarks. The form factor is computed at 
zero recoil using a new double ratio method that yields the form factor more directly than the 
previous Fermilab method. Other improvements over the previous calculation include the use of 
much lighter light quark masses, and the use of lattice (staggered) chiral perturbation theory in 
order to control the light quark discretization errors and chiral extrapolation. We obtain for the 
form factor, J-b^d*(X) = 0.921(13)(20), where the first error is statistical and the second is the 
sum of all systematic errors in quadrature. Applying a 0.7% electromagnetic correction and taking 
the latest PDG average for T B ->D*(l)\V cb \ leads to \V cb \ = (38.7 ± 0.9 e;cp ± 1.0^ eo ) x 1CT 3 . 

PACS numbers: 12.38.Gc, 13.25.Hw, 12.15.Hh 
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I. INTRODUCTION 



The Cabibbo-Kobayashi-Maskawa matrix element V c b plays an important role in the study 
of flavor physics Since \V c b\ is one of the fundamental parameters of the Standard 
Model, its value must be known precisely in order to search for new physics by looking 
for inconsistencies between Standard Model predictions and experimental measurements. 
For example, the Standard Model contribution to the kaon mixing parameter e# depends 
sensitively on \V c b\ (as the fourth power), and the present errors on this quantity contribute 
errors to the theoretical prediction of 6k that are around the same size as the errors due to 
Bx, the kaon bag parameter, which has been the focus of much recent work 
possible to obtain \V c b\ from both inclusive and exclusive semileptonic B decays, and both 
determinations are limited by theoretical uncertainties. The inclusive method fl, Q, fl 9. [fo| 
makes use of the heavy-quark expansion and perturbation theory. The method also requires 
non-perturbative input from experiment, which is obtained from the measured moments of 
the inclusive form factor B — > XJl>t as a function of the minimum electron momentum. The 
dominant uncertainties in this method are the truncation of the heavy quark expansion and 
perturbation theory ll|, ll2j. In order to be competitive with the inclusive determination 
of \V c b\ and thus serve as a cross-check, the exclusive method requires a reduction in the 
uncertainty of the B — > D* semileptonic form factor JF B ^ D », which has been calculated 
previously using lattice QCD in the quenched approximation [13]. 

Given the phenomenological importance of \V c b\, we have revisited the calculation of 
J-b^d* at zero recoil using the 2+1 flavor MILC ensembles with improved light staggered 
quarks [lj]] . The systematic error due to quenching is thus eliminated. The systematic 
error associated with the chiral extrapolation to physical light quark masses is also reduced 
significantly. Since staggered quarks are computationally less expensive than many other 
formulations, we are able to simulate at quite small quark masses; our lightest corresponds 
to a pion mass of roughly 240 MeV. Given the previous experience of the MILC Collaboration 
with chiral fits to light meson masses and decay constants [l^ ] , w e are in a regime where 
we expect rooted staggered chiral perturbation theory (rS^PT) |l7l. Ila . Il9l . [20j, [2l| to apply. 



We therefore use the rS^PT result for the B — > D* form factor [22J to perform the chiral 
extrapolation and to remove discretization effects particular to staggered quarks. In addition, 
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we introduce a set of ratios that allows us to disentangle light- and heavy-quark discretization 
effects, and we suggest a strategy for future improvement. Finally, we extract the B — > D* 
form factor using a different method from that originally proposed in Ref. This new 

method requires many fewer three-point correlation functions, and has allowed for a savings 
of roughly a factor of ten in computing resources, while at the same time simplifying the 
analysis. 

The differential rate for the semileptonic decay B — > D*tu^ is 

^ = ^s m U m B ~ m D *) 2 V^l G(w)\V cb \ 2 \F B ^D*(w)\ 2 , (1) 

where w = v' ■ v is the velocity transfer from the initial state to the final state, and 
Q{w)\J r B-*D* ( w )\ 2 contains a combination of four form factors that must be calculated non- 
perturbatively. At zero recoil G(l) = 1, and Tb^d*^) reduces to a single form factor, 
Given /i^fl), the measured decay rate determines \V c b\. 
The quantity h Al is a form factor of the axial vector current, 



{D*(v, e')\A^\B(v)) = W2m B 2m D . e^h Al {l) 



(2) 



where is the continuum axial- vector current and e' is the polarization vector of the D* . 
Heavy-quark s ym metry plays a useful role in constraining h Al (1), leading to the heavy-quark 
expansion [23|, \2A \ 

2£ A i P 



Va 



■v 



(2m c 



+ 



2m r 2mh 



(2m b y 



(3) 



up to order l/rr tn, a nd where r/A is a factor that matches heavy-quark effective theory 
(HQET) to QCD 



25 



261 ] . The €s are long distance matrix elements of the HQET. Heavy- 
quark symmetry forbids terms of order I/tuq at zero recoil 27|, and various methods have 
been used to compute the size of the 1/ttiq coefficients, including quenched lattice QCD 131 ] . 

The earlier work by Hashimoto et al. [3] used three double ratios in order to obtain 
separately each of the three l/m^ coefficients in Eq. ((3D. These three double ratios also 
determine three out of the four coefficients appearing at 1/ nig in the heavy-quark expansion. 



It was shown in Ref. [28j] that, for the Fermilab method matched to tree level in a s and to 
next-to-leading order in HQET, the leading discretization errors for the double ratios for this 
quantity are of order a s (A/2rriQ) 2 /g(amg) and (A/2mQ) 3 fi(amQ), where A is a QCD scale 
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stemming from the light degrees of freedom, such as that appearing in the HQET expansion 
for the heavy-light meson mass, rriM — m Q + A + ■ ■ • . The functions fi(arriQ) are coefficients 
depending on arriQ and a s , but not on A. When arriQ ~ 1, the fi(amq) are of order one; 
when arriQ <C 1, they go like a power of arriQ, such that the continuum limit is obtained. 
The powers of 2 are combinatoric factors. 

As discussed in Ref. 13], all uncertainties in the double ratios TZ used in that work 
scale as TZ — 1 rather than as TZ. Statistical errors in the numerator and denominator are 
highly correlated and largely cancel in these double ratios. Also, most of the normalization 
uncertainty in the lattice currents cancels, leaving a normalization factor close to one which 
can be computed reliably in perturbation theory. Finally, the quenching error, relevant to 



Ref. 



131 ] but not to the present unquenched calculation, scales as TZ — 1 rather than as TZ. 



This scaling of the error occurs because the double ratios constructed in Ref. 13| become 
the identity in the limit of equal bottom and charm quark masses. 

In the calculation reported here, the form factor /iAi(l) is computed more directly using 
only one double ratio, 

= { D*fr„m<Bfojy,c\D-) = a 

(D*\c~f 4 c\D*) (B\bj 4 b\B) V ; 

which is exact to all orders in the heavy-quark expansion in the continuum. 1 The lattice 
approximation to this ratio still has discretization errors that are suppressed by inverse 
powers of heavy-quark masses [a s (A/2mg) 2 and (A/2mq) 3 ], but which again vanish in the 
continuum limit. The errors in the ratio introduced in Eq. (jlj) do not scale rigorously as 
TZa± — 1 because TZa x is not one in the limit of equal bottom and charm quark masses. 
Nevertheless, this double ratio still retains the desirable features of the previous double 
ratios, i.e., large statistical error cancellations and the cancellation of most of the lattice 
current renormalization. Because the quenching error has been eliminated, the rigorous 
scaling of all the errors as TZ — 1, including the quenching error, is no longer crucial. The 
more direct method introduced here has the significant advantage that extracting coefficients 
from fits to HQET expressions as a function of heavy-quark masses is not necessary, and no 
error is introduced from truncating the heavy-quark expansion to a fixed order in l/rriq. In 



1 Note that the notation TZax stands for a different double ratio in Ref. 13 1 
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short, for an unquenched QCD calculation, the method using Eq. (j3J) gives a smaller total 
error than the method used in Ref. [jjj for a fixed amount of computer time . 

The currents of lattice gauge theory must be matched to the normalization of the contin- 



uum to obtain TZa 1 - The matching factors mostly cancel in the double ratio 29j, |30j, leaving 



= a/7^Ai = PyR~A\i where Ra 1 is the lattice double ratio and p, the ratio of matching 
factors, is very close to 1. (For the remainder of this paper we shall use the convention that 
a script letter corresponds to a continuum quantity, while a non-script letter corresponds 
to a lattice quantity.) This p factor has been calculated to one-loop order in perturbative 
QCD, and is found to contribute less than a 0.5% correction. We have exploited the p factors 
to implement a blind analysis. Two of us involved in the perturbative calculation applied 
a common multiplicative offset to the p factors needed to obtain /^(l) at different lattice 
spacings. This offset was not disclosed to the rest of us until the procedure for determining 
the systematic error budget for the rest of the analysis had been finalized. 

The unquenched MILC configurations generated with 2+1 flavors of improved staggered 
fermions make use of the fourth-root procedure for eliminating the unwanted four-fold degen- 
eracy of sta gger ed quarks. At non-zero lattic e sp acing, this procedure has small violations 



of unitarity 31 



32 



33 



34 



351 ] and locality 36J. Nevertheless, a careful treatment of the 



continuum limit, in which all assumptions are made explicit, argues that lattice QCD with 



rooted staggered quarks reproduces the desired local theory of QCD as a — > [37 



When 



39 



40 



411 ] for reviews) , this 



coupled with other analytical and numerical evidence (see Refs. 
gives us confidence that the rooting procedure is indeed correct in the continuum limit. 

The outline of the rest of this paper is as follows: Section II describes the details of the 
lattice simulation. Section III discusses the fits to the double ratios accounting for oscillating 
opposite-parity states. Section IV summarizes the lattice perturbation theory calculation of 
the p factor. Section V introduces the rooted staggered chiral perturbation theory formalism 
and expressions used in the chiral extrapolations. Section VI then discusses our treatment of 
the chiral extrapolation and introduces our approach for disentangling heavy and light-quark 
discretization effects. Section VII provides a detailed discussion of our systematic errors, and 
we conclude in Section VIII. 
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II. LATTICE CALCULATION 



The lattice calculation was done on the MILC ensembles at three lattice spacings with 
a ~ 0.15, 0.125, and 0.09 fm; these ensembles have an 0(a 2 ) Symanzik improved gauge 



action and 2+1 flavors of "AsqTad" improved staggered sea quarks [42|, |43j, |4j, |45j, |46|, 1471 ]. 
The parameters for the MILC lattices used in this calculation are shown in Table [H We 
have several light masses at both full QCD and partially-quenched points (m va \ ^ m sea ), 
and our light quark masses range between m s /10 and m s /2. Table HT1 shows the valence 
masses computed on each ensemble. In this work we follow the notation where m s is 
the physical strange quark mass, m is the average u-d quark mass, and m', m! s indicate 
the nominal values used in simulations. In practice, the MILC ensembles choose m' s within 
10-30% of m s and a range of fh' to enable a chiral extrapolation. 

The heavy quarks are computed using the Sheikholeslami-Wohlert (SW) "clover" action 



481 ] with the Fermilab interpretation via HQET 49] . The SW action includes a dimension- 
five interaction with a coupling csw that has been adjusted to the value Mq 3 suggested by 



tadpole- improved, tree-level perturbation theory [50j. The value of Uq is calculated either 
from the plaquette (a ~ 0.15 fm and a ~ 0.09 fm), or from the Landau link (a ~ 0.12 fm). 
The adjustment of csw is needed to normalize the heavy quark's chromomagnetic moment 
correctly 49]. 

The tadpole-improved bare quark mass for SW quarks is given by 

u \2k 2K clit J 

where tuning the parameter k to the critical quark hopping parameter K crit would lead to 
a massless pion. The spin averaged B s and D s kinetic masses are computed on a subset 
of the ensembles in order to tune the bare k values for bottom and charm (and hence the 
corresponding bare quark masses) to their physical values. These tuned values were then 
used in the B — > D*lv form-factor production run. 

The relative lattice scale is determined by calculating r\ja on each ensemble, where r\ is 



related to the force between static quarks by r\F{r\) = 1.0 5lJ, |52j. To avoid introducing 
implicit dependence on m', m' s via ri(m' ,m' s , g 2 ) (where, as above, primes denote simulation 
masses), we interpolate in m' s and extrapolate in m' to obtain ri(m, m s , g 2 )/a at the physical 
masses. We then convert from lattice units to r\ units with ri(fh,m s , g 2 )/a. Below we shall 



TABLE I: Parameters of the simulations. The columns from left to right are the approximate lattice 
spacing in fm, the sea quark masses arh'/am' s , the linear spatial dimension of the lattice ensemble 
in fm, the dimensionless factor m^L (m^ corresponds to the taste-pseudoscalar pion composed of 
light sea quarks), the gauge coupling, the dimensions of the lattice in lattice units, the number of 
configurations used for this analysis, the bare hopping parameter used for the bottom quark, the 
bare hopping parameter used for the charm quark, and the clover term csw used for both bottom 
and charm quarks. 



a(fm) 


afh'/am's 


L(fm) 




10/ff 2 


Volume # Confi| 


JS K b 


K c 


csw 


0.15 


0.0194/0.0484 


2.4 


5.5 


6.586 


16 3 x 48 


628 


0.076 


0.122 


1.5673 


0.15 


0.0097/0.0484 


2.4 


3.9 


6.572 


16 3 x 48 


628 


0.076 


0.122 


1.5673 


0.12 


0.02/0.05 


2.4 


6.2 


6.79 


20 3 x 64 


460 


0.086 


0.122 


1.72 


0.12 


0.01/0.05 


2.4 


4.5 


6.76 


20 3 x 64 


592 


0.086 


0.122 


1.72 


0.12 


0.007/0.05 


2.4 


3.8 


6.76 


20 3 x 64 


836 


0.086 


0.122 


1.72 


0.12 


0.005/0.05 


2.9 


3.8 


6.76 


24 3 x 64 


528 


0.086 


0.122 


1.72 


0.09 


0.0124/0.031 


2.4 


5.8 


7.11 


28 3 x 96 


516 


0.0923 


0.127 


1.476 


0.09 


0.0062/0.031 


2.4 


4.1 


7.09 


28 3 x 96 


556 


0.0923 


0.127 


1.476 


0.09 


0.0031/0.031 


3.4 


4.2 


7.08 


40 3 x 96 


504 


0.0923 


0.127 


1.476 



call this procedure the mass-independent determination of r\. 

In order to fix the absolute lattice scale, one must compute a physical quantity that can 



be compared directly to experiment; we use the T 2S-1S splitting [53j and the most recent 
MILC determination of f n 54|. The difference between these determinations results in a 
systematic error that turns out to be much smaller than our other systematics. When the 



T scale determination is combined with t 
ue rf ys = 0.318(7) fm 



quark masses, a va 
0.3108(15)(^)fm 



ie continuum extrapolated r\ value at physical 
55] is obtained. The f n determination is r^ hys = 



54]. Given r^ hys , it is then straightforward to convert quantities measured 



in T\ units to physical units. 

The dependence on the lattice spacing a is mild in this analysis. Since a only enters the 
calculation through the adjustment of the heavy and light quark masses, the dependence of 
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TABLE II: Valence masses used in the simulations. The columns from left to right are the approx- 
imate lattice spacing in fm, the sea quark masses am! jam! s identifying the gauge ensemble, and the 
valence masses computed on that ensemble. 



o(fm) 


Clffl' 1 CLTTl 


a,Tn x 


» 0.15 


0.0194/0.0484 


0.0194 


ps 0.15 


0.0097/0.0484 


0.0097, 0.0194 


ps 0.12 


0.02/0.05 


0.02 


ps 0.12 


0.01/0.05 


0.01, 0.02 


ps 0.12 


0.007/0.05 


0.007, 0.02 


ps 0.12 


0.005/0.05 


0.005, 0.02 


ps 0.09 


0.0124/0.031 


0.0124 


ps 0.09 


0.0062/0.031 


0.0062, 0.0124 


ps 0.09 


0.0031/0.031 


0.0031, 0.0124 



ft-Ai(l) on a is small. Staggered chiral perturbation theory indicates that the a dependence 
coming from staggered quark discretization effects is small 22[], and this is consistent with 
the simulation data. 

In this work, we construct lattice currents as in Ref. 



49] 



J, 



hh' 



z%z" h 'y h r^ h ,, 



(6) 



where is either the vector (z7 M ) or axial-vector (27^75) current. The rotated field is 
defined by 



*/, = (! + arf l7 • D lat )^ 



(7) 



where iph is the (heavy) lattice quark field in the SW action. Di at is the symmetric, nearest- 
neighbor, covariant difference operator; the tree-level improvement coefficient is 



1 



uq \2 + m^a 2(l + moa) 

In Eq. ([6]) we choose to normalize the current by the factors of Zy^ [h = c,b) since even 
for massive quarks they are easy to compute non-perturbatively. The continuum current is 
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related to the lattice current by 



J, 



hh' 



up to discretization effects, where 



n T hh ' 



ybc *ycb 
Jr Jr 
ycc ybb 



(9) 



(10) 



bb ryoc 



and the matching factors Zj^' 7 s are defined in Ref. 30]. Note that the factor yZ^Z 
multiplying the lattice current in Eq. ([6]) cancels in the double ratio by design, leaving only 
the p factor, which is close to one and can be computed reliably using perturbation theory. 
The perturbative calculation of pj r is described in more detail in Section IV. 

Interpolating operators are constructed from four-component heavy quarks and staggered 
quarks as follows. Let 



ifj b (x)~f 5 n(x)x(x), 



(11) 

(12) 



where x is the one-component field in the staggered- quark action, and 



xi/a x 2 /a x 3 /a x^/a 

VL{x) = 7! 7 2 7 3 7 4 . 



(13) 



The left (right) index of (Q) can be left as a free taste index 41(] or x can be promoted 
to a four-component naive-quark field to contract all indices 56J]. The resulting correlation 
functions are the same if the initial and final taste indices are set equal and then summed. 
The same kinds of operators have been used in previous calculations 



57 



58 



59J. 



Lattice matrix elements are obtained from three-point correlation functions. The three- 
point correlation functions needed for the B — > D* transition at zero-recoil are 



C B ^ D * (U, t„t f ) = ( x ' */)*c7;75* 6 (y, t.)Ol(0, U)\0), 

C B ^ B (U,t s ,t f ) = ^<0|O fl (x,t / )W 6T ^ 6 (y,t a )O^(0,t i )|0), 
C D *^ D *{ti,t s ,t f ) = ^(0|C? D .(x, t f )q cl ^ c (y,t s )O^(0,U)\0). 



(14) 
(15) 
(16) 



x,y 
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In C B ^ D * the polarization of the D* lies along spatial direction j. If the source-sink separa- 
tion is large enough then we can arrange for both t s — ti and tf — t s to be large so that the 
lowest-lying state dominates. Then 

y/2m D *y/2m B 

where and m^. are the masses of the B and D* mesons and Zh = \ (0\Oh\H)\ 2 . 

In practice, the meson source and sink are held at fixed tj = and tf = T, while the 
operator time t s = t is varied over all times in between. Using the correlators defined in 
Eqs. ( 1141116]) we form the double ratio 

C B - D '(0,t,T)C D '- B (0,t,T) 
Al{) C D ^ D *(0,t,T)C B ^ B (0,t,T)' { ] 



All convention-dependent normalization factors, including the factors of \J Zjijlrxiu-, cancel 
in the double ratio. In the window of time separations where the ground state dominates, a 
plateau should be visible, and the lattice ratio is simply related to the continuum ratio *R-Ax 
by a renormalization factor 



pa 1 \/Rm = VnZ = h Al (i), (19) 

with as in Eq. (TTOl) . The right-hand side of Eq. fTTTT) is the first term in a series, with 
additional terms for each radial excitation, including opposite-parity states that arise with 
staggered quarks. Eliminating the opposite-parity states requires some care, and this is 
discussed in detail in the next section. In order to isolate the lowest-lying states we have 
chosen creation and annihilation operators, 0^ B and Op*, that have a large overlap with 
the desired state. This was done by smearing the heavy quark and anti-quark propagator 
sources with IS Coulomb-gauge wave-functions. 

III. FITTING AND OPPOSITE-PARITY STATES 

Extracting correlation functions of operators with staggered quarks presents an extra 
complication because the contributions of opposite-parity states introduce oscillations in 
time into the correlator fits 



561 ] . Three-point functions obey the functional form 

-ir(-l)^-'A ek e— 

k=0 1=0 

11 



C x ^ Y (0,t,T) = ^^-l^-l^-'^e^V^-'). (20) 



For odd k and I the excited state contributions change sign as the position of the operator 
varies by one time slice. Although they are exponentially suppressed, the parity partners 
of the heavy-light mesons are not that much heavier than the ground states in which we 
are interested, so the oscillations can be significant at the source-sink separations typical of 
our calculations. These separations cannot be too large because of the rapid decrease of the 
signal due to the presence of the heavy quark. 

Although one can fit a given three-point correlator to Eq. ()2Up . in the calculation of 
/^(l) we use double ratios in which numerator and denominator are so similar that most of 
the fitting systematics cancel, and it is convenient to preserve this simplifying feature. We 
do this by forming a suitable average over correlator ratios with different (even and odd) 
source-sink separations. It turns out that the amplitudes of the oscillating states in B — > D* 
correlation functions are much smaller than they are in many other heavy-light transitions 
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and that the oscillating states in B — > D* are barely visible at the present level 
of statistics. Even so, we introduce an average that reduces them still further, to the point 
where they are negligible. 

Although we shall take the average of the double ratio, let us first examine the average of 
an individual three-point function. Expanding Eq. ( 1201) so that it includes the ground state 
and the first oscillating state, we have 

+ (-l)*A^ y e" m ^ t "' mr(T "* ) + (-l) T Af 1 ^ y e~ m * i ~ m ' y(T ~ t) + ... 

= A x-+Y e - mx t- mY {T-t) ^ + c *^(o, t, T) + ...], (21) 

where in the last line we have pulled out the ground state amplitude and exponential de- 
pendence. The function c x_>y (0,t,T) is defined 

C X ^ y (0,t,T) EE if^(-lf-V Am ^ T -') + if^(-l)V Am *' 
^00 ^00 

+ ^_(_l)^-^-Am y (T-t) ) (22) 

where Amx,Y = m 'xY ~ m x,Y is the splitting between the lowest-lying desired-parity state 
and the lowest-lying wrong-parity state. Note that the first two terms produce oscillations 
as the position of the operator is varied over the time extent of the lattice. The third 
term, however, changes sign only when the total source-sink separation is varied. It is this 

12 



term that our average is designed to suppress, since it will not be as clearly visible in the t 
dependence of the lattice data as those that oscillate in t. 
We define the average to be 



C (0,t,T) 



±C x ^ Y (0,t,T) + ~C x ^ Y (0,t,T + l) 



C x ^ Y (0,t+l,T + l) 



(23) 



~^ y (0,t,T) from Eq. (121]) into this definition gives 
C A ^ r (0, t, T) = A$ J -* Y e- mxt - mY V- t) [l + c x ~* y (0, t, T) + ...], (24) 



Substituting the expression for C 

-^X^Y 



where the function c 



X^Y : 



-X^Y 



(0,t,T) 



IS 

AX^Y 

AX->Y\ 
^00 

A X ~* Y 
+ 10 | 

^oo 

X^Y 



_-|\T-t e -Amy(T-i) 



1 1, 

2 + 4^ 



-Amy > 



■i)V 



4 A 



^oo 



^^Tg-Amxt-AmytT-i) 



1 _ l( e - A »"y _|_ g-Am x N 



(25) 



Note that Eq. (125]) has the same exponential time dependence as Eq. (1221) . but with the size 
of the amplitudes reduced by the factors in square brackets. Thus, the average is equivalent 
to a smearing that reduces the oscillating state amplitudes. It is possible to compute the 
Arrix precisely from fits to two-point correlators. We find values between about 0.2 and 0.4 
in lattice units. Given these values, the first two factors in brackets reduce their respective 
amplitudes by approximately a factor of 2, and the targeted, non-oscillating term is reduced 
by a factor of ~ 6-10. 

Specializing to the B — > D* case, consider the double ratio 



^00 ^00 



— c 



D*—>D* i 



'(0,t,T) - c M (0, t,T) + ...], (26) 

where we have again factored out the ground state contribution. Equation (126]) follows from 
Eq. (j!8p treating the c's as small. Note that the c's are expected to be similar in numerator 
and denominator, and to the extent that they are the same they will cancel in this expression. 
Applying the average in Eq. (T23]) directly to the double ratio, 

R(0,t,T) = ±R(0,t,T) + ±R(0,t,T+l) 

(27) 



+ - J R(0,t+l,T + l), 



13 




FIG. 1: Double ratio Ra 1 on the am' = 0.0124 fine (a = 0.09 fm) ensemble. The source was fixed 
to time slice 0, and the operator position was varied as a function of time. Two different sink points 
were used with even and odd time separations between source and sink in order to study the effect 
of non-oscillating contributions from wrong parity states. 



we get 



RM(0,t,T) = SZ A Xl [l + c»- D \0,t,T) + c D -»*(0,t,T) 

^inn Jinn 



"•oo 1 ^oo 



(0, t,T) (0,t,T) + ...], (28) 

where each of the oscillating state terms in the individual three-point functions is suppressed 
according to Eq. (125]) . 

Although Arris and Amp* can be obtained from fits to the two point correlators, the 
oscillating state amplitudes appearing in the three-point correlators must be determined 
directly from the three-point correlator data. Figure Q] shows the double ratio Ra 1 used to 
obtain /^(l). The source is at time slice 0, the sink is at T, and the operator position is 
varied along t. Two different source-sink separations were generated that differed by a single 
time unit at the sink (T = 17, 18). The average of these two correlators was taken according 
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FIG. 2: Averaged double ratio, Ra 1} of Eq. (|27p on the am' = 0.0124 fine (a = 0.09 fm) ensemble. 
The plateau fit is shown with la error band. 

to Eq. fl27|) . and this average was fit (including the full covariance matrix) to a constant, as 
shown in Fig. El There is no detectable oscillation even before the average is taken, as can 
be seen in Fig. [TJ according to Eq. (pZol) the oscillating contributions are reduced even further 
in the average so that their systematic errors can be safely neglected. 



IV. PERTURBATION THEORY 



Lattice perturbation theory is needed in order to calculate the short-distance coefficient 
Pa 1 defined in Eq. f fTUj) . Although naive lattice perturbation theory appears to converge 
slowly, the two main causes have been identified jsoj: the bare gauge coupling is a poor 
expansion parameter, and coefficients are large when tadpole diagrams occur. If a renormal- 
ized coupling is used as an expansion parameter, and one computes only those quantities 
for which the tadpole diagrams largely cancel, then lattice perturbation theory seems to 
converge as well as perturbation theory in continuum QCD. 
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Only the vertex correction contributes to the p factor, as the wave-function renormaliza- 
tion (including all tadpoles) cancels by construction. Even the vertex correction partially 
cancels, and the one-loop coefficient is found to be small. The perturbative corrections to 
the p factor can be written as 



hh' 
PJt 



l + a y (g*)47rp^ [1] + - 



yhh yh'h' 



(29) 



where p h ]^ ^ is the coefficient of the one-loop correction, and the coupling ay is the renor- 



malized strong coupling constant in the V-scheme 5(3, |62j], which is based on the static-quark 



potential. The coupling is determined following the procedure of Ref. 



63|. The scale q* of 



the running coupling ay(q*) should be chosen to be the typical momentum of a gluon in 



the loop. A prescription 



Mackenzie (BLM) 
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or calculating this scale was introduced by Brodsky, Lepage, and 
62] . They define q* by 

Jd 4 qf(q) ln(g 2 ) 



ln(g 



*2 > 



(30) 



fd 4 qf(q) 

where f(q) is the one- loop integrand and the numerator is the first log moment. This pre- 
scription was extended by Hornbostel, Lepage, and Morningstar (HLM) 64| to cases where 
the one-loop contribution is anomalously small leading to a break down of Eq. (I3"0l . The 
HLM prescription for q* takes into account two-loop contributions to the gluon propagator 
via the inclusion of second log moments. Since we do encounter anomalously small one-loop 
corrections in p^, the HLM prescription was used to determine q*. Results for q* HLM and 
Pa x needed for this calculation are given in Table IIHI The p factor varies somewhat as a 
function of lattice spacing, and is even slightly different from ensemble to ensemble at the 
same nominal lattice spacing, due to the slightly different (5 values used to generate the 
gauge fields. 



The calculation of pA x is described in Refs. 65|, |66j . It uses automated perturbation theory 



techniques to generate the Feynman rules and VEGAS [67( for the numerical integration of 

the loop integrals. As a check, it was verified that this calculation reproduces known results 

□ 

for the heavy-heavy currents with the Wilson plaquette action [29J and for the V4 current in 
the massless limit with the Symanzik improved gauge action. 

As mentioned in the introduction, we have exploited the p factor to implement a blind 
analysis. Two of us applied a multiplicative offset close to 1 to the p factor, generated with a 
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TABLE III: Computed values of pa 1 in the HLM prescription [64| . The first three columns label 
each ensemble with the approximate lattice spacing in fm, the light sea quark mass am', and the 
strange quark mass am' s . The fourth column is aq* HLM , where the error is calculated using the 
statistical error from VEGAS for the Oth, 1st, and 2nd moments of the one-loop integrals. The 
fifth column is pa 1 on that ensemble, and the errors are the statistical errors from the VEGAS 
evaluation, including the one-loop coefficients and q* HLM . 



(fm] 


am! 


am' s 


a 1*HLM 




0.15 


0.0194 


0.0484 


2.03(10) 


0.9966(2) 


0.15 


0.0097 


0.0484 


2.03(10) 


0.9966(2) 


0.12 


0.02 


0.05 


1.96(10) 


0.9964(2) 


0.12 


0.01 


0.05 


1.96(10) 


0.9964(2) 


0.12 


0.007 


0.05 


1.96(10) 


0.9964(2) 


0.12 


0.005 


0.05 


1.96(10) 


0.9964(2) 


0.09 


0.0124 


0.031 


2.98(14) 


1.00298(9) 


0.09 


0.0062 


0.031 


2.98(14) 


1.00300(9) 


0.09 


0.0031 


0.031 


2.98(14) 


1.00301(9) 



random key. The offset was not unlocked until the procedure for determining the systematic 
errors in the rest of the analysis had been finalized. 

V. STAGGERED CHIRAL PERTURBATION THEORY 

The simulation masses m' val and m' sea (for valence and sea) are all larger than the phys- 
ical m. A controlled chiral extrapolation can be guided by an appropriate chiral effective 
theory that includes the effect of staggered- quark discretization errors. Rooted staggered 



chira 



perturbation theory (rS^PT), which has been formulated for heavy-light quantities in 



Ref. 68J, is such a theory. In rS^PT, a replica method is used to take into account the effect 



of rooting; this procedure has been justified in Refs. [33|, |6 

Because of taste-symmetry breaking, the staggered theory has 16 light pseudoscalar 
mesons instead of 1. The tree- level relation for the masses of light staggered mesons in 
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the chiral theory is [l?], 



m xy,E = Vo(m x + m y ) + a An , (31) 



where and m y are staggered quark masses, /xq is the continuum low-energy constant, and 
a 2 A= are the splittings of the 16 pions of taste S. For staggered quarks there exists a residual 
SO (4) taste symmetry broken at O (a 2 ), such that there is some degeneracy among the 16 



pions [17|, and the taste index H runs over the multiplets P, A,T,V, I with degeneracies 1, 
4, 6, 4, 1. The splitting a 2 Ap vanishes because there is an exact (non-singlet) lattice axial 
symmetry. 

Schematically, the next-to-leading order (NLO) result for the relevant form factor is 

h%°0)/VA = 1 + X A {K X ) + Jj^ x log Sl „ loop (A x ) (32) 

where X^(A X ) is a low energy constant of the chiral effective theory, and is therefore indepen- 
dent of light quark mass and cancels the chiral scale dependence A x of the chiral logarithms. 
By heavy-quark symmetry, Xa{A- x ) is proportional to 1/m 2 in the heavy-quark expansion. 
The term i]a is a factor which matches heavy-quark effective theory to QCD, and contains 
perturbative-QCD logarithmic dependence on the heavy-quark masses; it is independent of 
light quark mass. The term proportional to g\j)* v is short-hand for the one-loop staggered 
chiral logarithms, and is given in the appendix for ease of reference. The rooted staggered ex- 



pression was derived in Ref. [22J. The one- loop staggered logarithms depend on both valence 
and sea quark masses, and include taste-breaking effects coming from the light quark sector. 
This expression also contains explicit dependence on the lattice spacing a, and requires as 
inputs the papers of the staggered chiral W * A , m addition to the staggered 
taste splittings &p,a,t,v,i These parameters can be obtained from chiral fits to the 

light pseudoscalar meson sector and are held fixed in the chiral extrapolation of /^(l). The 
continuum low energy constant gDD*n appears, and below we take a generous range inspired 
by a combined fit to many different experimental inputs, including a leading-order analysis 
of the D* width. The D*-D splitting is well determined from experiment. The only 
other parameter that appears at NLO is the constant Xa(A), and this is determined by our 
lattice data for /^(l). 

Although the lattice data are well described by the NLO formula, it is useful to go 
beyond NLO and to include the next-to-next-to-leading-order (NNLO) analytic terms as a 
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way to estimate systematic errors. We do not include the NNLO logarithms because they 
are unknown and would require a two-loop calculation. The expression including analytic 
terms through NNLO is 

h^ LO (l)/ VA = 1 + NLO + Cl m 2 Xp + c 2 (2m 2 Up + m 2 Sp ) + c 3 a 2 , (33) 

where the subscript P on the meson masses indicates the taste pseudo-scalar mass. We use 
the notation from the rS^PT literature that mx s is a taste H meson made of two valence x 
quarks, mt/ a is a taste H meson made of two light sea quarks, and ms s is a taste S meson 
made of two strange sea quarks. By heavy-quark symmetry, the q are suppressed by a factor 
of 1/m 2 . Since the only free parameter through NLO is an overall constant, we include the 
NNLO analytic terms in the fit used for our central value. This leads to a larger statistical 
error and is more conservative. 



VI. TREATMENT OF CHIRAL EXTRAPOLATION 



In this section, we discuss the approach we have developed to disentangle the heavy- and 
light-quark discretization effects and to perform the chiral and continuum extrapolations. 
In the Fermilab method, heavy-quark discretization errors can be estimated by comparing 



the heavy-quark expansions for lattice gauge theory and continuum QCD [28|, |29j, |30|, 1 701 ]. 
The dependence on a is not simply a power series (unless ma <C 1), so power-counting 
estimates in HQET are used. On the other hand, some of the light quark discretization 
effects are constrained by rS^PT. The heavy-quark errors are asymptotically constrained 
by the Symanzik low-energy Lagrangian when m^a <C 1 and by heavy-quark symmetry 
even when m^a is close to 1. In the region in between, the errors smoothly interpolate the 



asymptotic behavior 49 



70[ . The errors in the SW action used for the heavy quarks decrease 



with lattice spacing as a s a in the m^a <C 1 region, as compared to the light quark (improved 
staggered) discretization errors, which decrease much faster, as a s a 2 . 

The first step of the method is to normalize the numerical data for h Al (l) to a fiducial 
point by forming the ratio 

h Al (m x ,m',m' s ,a) 



n M (m x ,m ,m.,a) 



/i Al (mfi d ,m fid ,mf d ,a) 



(34) 
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where m fid is a fiducial mass, m x is the light (spectator) valence quark, fh' is the isospin 
averaged light sea quark on a particular ensemble, and m' s is the strange sea quark on that 
ensemble. (Note that the factor of t)a in Eqs. (1321) and (1331) cancels in the ratio.) The principle 
advantage of this ratio is that heavy quark discretization effects largely cancel, since the heavy 
quarks are the same in numerator and denominator. This allows us to disentangle the heavy- 
quark discretization effects from those of the light quark sector coming from staggered chiral 
logarithms, thus isolating the (taste-violating) discretization effects specific to the staggered 
light quarks. These light quark discretization effects can appear in non-analytic terms in 
rS^PT and are due to violations of taste-symmetry. They can be removed to a given order 
in rS^PT (we work to NLO) in fits to the numerical data at multiple lattice spacings using 
the explicit rS^PT formula of Eq. ( l33l) . since this formula includes the staggered lattice 
artifacts. The continuum limit of the ratio TZm can be obtained using our fitted values for 
parameters in rS^PT and taking a — > in the rS^PT expression for TZm- We do not need a 
more explicit ansatz for the functional form of the heavy-quark discretization effects, since 
they largely cancel in the ratio. 

Normalizing the continuum extrapolated ratio TZ^ by 1ia 1 at the fiducial point on a very 
fine fiducial lattice where the heavy-quark discretization effects are small gives a value close 
to the physical continuum result, 

hA 1 {rh,m,m s ,0) « ^(mf ,m fid ,mf ,a fid ) x 1Z&<i(m, fh, m s , 0), (35) 

where the relation becomes exact as a fid — > 0. Note that the requirement that the heavy- 
quark discretization effects must be small enforces the condition that the improved staggered 
light-quark discretization effects be even smaller (and likely negligible) because the staggered 
discretization effects decrease much faster with lattice spacing. The fiducial masses m^ d , m fid , 
and rrig should be chosen large enough that it would be feasible to simulate this mass point 
on a very fine lattice (since the cost rises significantly as the mass of the light sea quarks is 
decreased), thus normalizing the lattice data to a point where the heavy-quark discretization 
effects are small. The fiducial masses should not be chosen so large, however, that rS%PT 
would not be a reliable guide in performing the continuum and chiral extrapolation of 71m- 
This method can be considered the crudest form of step-scaling, but it does illustrate that 
one does not need lattices which are simultaneously fine enough for b quarks and large enough 
for light quarks in order to simulate, with high precision, quantities that involve both. In 
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practice, we find m^ d = m fid « 0.4m s and m^ d m m s are reasonable values for the fiducial 
masses. The fiducial lattice spacing should be chosen as fine as is practical; a succession of 
progressively finer fiducial lattices would be desirable for verifying that the a dependence 
is of the expected size. In this work we take our finest lattice (0.09 fm) as the fiducial 
lattice, but we apply Eq. (J35l) with the coarser lattices taken as fiducial lattices in order to 
estimate discretization errors. We note that the method presented above can be applied to 
all calculations involving the Fermilab treatment of heavy-quarks and staggered light quarks, 
not only the B — > D*£u form factor h Al . It may also be desirable to compute quantities 
at the fiducial point (or a succession of such points) using an even further improved action 
for the heavy quarks. Once the fiducial lattice spacing is of the order 0.03-0.01 fm, even 
the bottom quark may be treated as a "light" quark with the highly improved staggered 



action (HISQ) [71[ or with chiral fermions, for which mass dependent discretization effects 
are small. Conserved currents could then be used for many simple heavy-light quantities, 
removing the need for a perturbative renormalization. 

For the chiral extrapolation of h Al we find it useful to form two additional ratios, 

ftsea(m ,m s ,a) = ^— , 36 

h Al {m% a , m ba , m* d , a) 

v ( -/ / \ _ h Al (m x ,m',m' a ,a) 

U Y ^{m x ,m ,m s ,a) = - — - — — - — r, (37) 

h Al {m ba ,m f ,m' s ,a) 



whose product is clearly 71m, Eq. (1311) . 7^ sea and 7£ va i separate the sea and valence quark 
mass dependence, which makes it easier to assess systematic errors. The values of h Al that 
enter Eqs. (1361) and (1371) are obtained from 



h Al =pyRA 1 , (38) 

where R Al is the average of double ratios defined in Eqs. (1281) . The ratios in Eqs. (136|) and (1371) 
are now quadruple ratios, where the excited state contamination is further suppressed over 
that of the double ratio. Performing the chiral extrapolation, taking the continuum limit of 
the two ratios, and multiplying them together we recover 7^fid(^ ; m, m s, 0) by construction. 
Thus, we can rewrite Eq. (1351) as 

< ys « h Al {mf, m fid , mf, a M ) x [ft sca (m, m 9 , 0) x ^ val (m, m, m s , 0)], (39) 

where, again, the relation becomes exact as a fid — > 0. 
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TABLE IV: Fiducial masses used at the three different lattice spacings. The first four columns are 
the approximate lattice spacing in fm, the fiducial valence quark mass, the fiducial light sea quark 
mass, and the fiducial strange quark mass. The fifth and sixth columns are the values of \l Ra ± 
and , respectively, computed at that fiducial point. 

lattice spacing (fm) amf am fid amf h% 

0.15 0.0194 0.0194 0.0484 0.9211(73) 0.9180(73) 

0.12 0.02 0.02 0.05 0.9112(73) 0.9079(73) 

0.09 0.0124 0.0124 0.031 0.9210(85) 0.9237(85) 



To the extent that the extrapolation in sea quark masses is mild, the ratio TZ sca should be 
close to one, since the valence light mass is the same in both numerator and denominator. 
lZ va \ contains less trivial chiral behavior. However, since the numerator and denominator are 
computed on the same ensemble (with different valence masses), they are correlated, and 
statistical errors tend to cancel in TZ va \. The ratio TZ sea , has small statistical errors because the 
valence mass m^ d in that ratio is relatively heavy. Of course, the heavy-quark discretization 
errors are significantly suppressed in both ratios, isolating the light quark mass dependence 
and staggered discretization effects. A direct chiral fit to the numerical data (not involving 
the ratios introduced here) would require a more explicit ansatz for the treatment of the 
heavy quark discretization effects than is needed in the ratio fits 2 . Note that in the ratios 
the fiducial point need not be tuned to the same mass at every lattice spacing; differences 
can be accounted for in the fit itself. The fiducial points used at different lattice spacings 
are m^ d = m fid = 0Am' s and mf d = m' s . The explicit values are given in Table IIV] along 
with the calculated values of \J Ra 1 and /i^ d at that fiducial point. 

The constant term A^(A X ) in Eq. (1521 cancels in the ratios 7Z sea and 7£ va i, so the behavior 
of these ratios is completely predicted through NLO in the chiral expansion. We find good 
agreement between the predicted form and the numerical data. However, given that our 
fiducial spectator quark mass is rather large (around 0.4m s ), we include the NNLO analytic 

2 A direct (correlated) chiral fit would still, however, reflect the correlations which cause cancellations in 
the statistical errors in the ratios. 
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terms in the ratio fits in order to estimate systematic errors associated with the chiral 
expansion. There are only two new continuum low energy constants introduced at this 
higher order, and the ratios TZ sea . and TZ va \ determine one each. There is also an analytic 
term proportional to a 2 appearing at this order, but it cancels in each of the TZ sca and TZ va \ 
ratios. 

In future calculations, it would be feasible to use a much finer lattice spacing for the 
fiducial point, thereby further reducing heavy-quark discretization errors. For now, however, 
we use hA 1 (fn^ d , m fid , mf d , 0.09 fm), with the fiducial masses in Table HV] in Eq. (|39|) . As 
a way to estimate discretization errors we use our results for h^, at the two coarser lattice 
spacings in Eq. (|39|) also. 

At the lattice spacings used in this work the light-quark discretization effects may still 
be non-negligible compared to heavy-quark discretization effects. With rS^PT it is possible 
to remove from the discretization effects associated with staggered chiral logarithms, 
although purely analytic discretization errors remain. Removing this subset of staggered 
effects leads to a value for the fiducial form-factor which we call the "taste-violations-out" 
value. Not removing them leads to the "taste-violations-in" value. The difference turns out 
to be negligible, less than 0.1% on our coarsest ensemble and less than 0.01% on the fine 
ensemble. Thus, the discretization effects in our lattice data coming from taste-violations 
in the staggered chiral logarithms are extremely small at the fiducial point mass, and we 
neglect this difference in the analysis. 

Figure [3] shows the plateau fit to the ratio TZ va \ on the fine ensemble with (am',am' s ) = 
(0.0062,0.031). The valence mass in the numerator is the full QCD value of am' x = 0.0062, 
while the fiducial valence mass in the denominator is am^ d = 0.0124. Both numerator 
and denominator are computed on the same ensemble, so they have the same sea quark 
masses, and correlated statistical errors largely cancel in the ratio, as expected. Excited- 
state contamination is also reduced. Computed values for TZ sca on all of our ensembles are 
given in Table [V] and the computed values for 7£ val are given in Table [VTl 
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1.06 



1.04 



1.02 



0.98 



0.96 



0.94 




FIG. 3: 7^ va i on the am! = 0.0062 fine ensemble. The valence mass in the numerator is the full 
QCD value of am' x = 0.0062 while the fiducial valence mass in the denominator is am^ = 0.0124. 
The fit to a constant has a x 2 /d.o.f = 0.20. 



TABLE V: Computed values of TZ se3 .. The first three columns are the arguments of TZ sea as defined 
in Eq. ()36[) : they are the light sea quark mass m', the strange quark mass m' s , and the approximate 
lattice spacing in fm. The fourth column is 7£ sea . 



am! 


am' s 


a(fm) 




0.0097 


0.0484 


0.15 


1.009(12) 


0.01 


0.05 


0.12 


1.0070(98) 


0.007 


0.05 


0.12 


1.0027(91) 


0.005 


0.05 


0.12 


1.014(10) 


0.0062 


0.031 


0.09 


1.000(12) 


0.0031 


0.031 


0.09 


0.996(10) 
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TABLE VI: Computed values of TZ VSL \. The first four columns are the arguments of lZ va \ as defined 
in Eq. (|37p ; they are the light valence quark mass m x , the light sea quark mass in', the strange 
quark mass m' s , and the approximate lattice spacing in fm. The fifth column is TZ VSb \. 



am x 


am! 


am' 

S 


a(fm) 




0.0097 


0.0097 


0.0484 


0.15 


1.0056(65) 


0.01 


0.01 


0.05 


0.12 


0.9994(41) 


0.007 


0.007 


0.05 


0.12 


0.9900(57) 


0.005 


0.005 


0.05 


0.12 


1.0081(90) 


0.0062 


0.0062 


0.031 


0.09 


1.0005(50) 


0.0031 


0.0031 


0.031 


0.09 


1.0043(62) 



VII. SYSTEMATIC ERRORS 

In the following subsections, we examine the uncertainties in our calculation due to fit- 
ting and excited states, the heavy-quark mass dependence, the chiral extrapolation of the 
light spectator quark mass, discretization errors, and perturbation theory. As mentioned in 
Section II, statistical uncertainties are computed with a single elimination jackknife and the 
full covariance matrix. 

A. Fitting and excited states 

We have examined plateau fits to the time dependence of the double and quadruple ratios 
introduced in Sections I and V. The \ 2 m our fits is defined with the full covariance matrix. 
The fits to the ratios were done under a single elimination jackknife, after blocking the 
numerical data by 8 on the fine lattices and by 4 on the coarse and coarser lattices. The 
blocking procedure averages 4 (or 8) successive configurations before performing the single 
elimination jackknife. These values for the block size were chosen such that the statistical 
error on the double ratio fit did not increase when a larger block size was used. Statistical 
errors were determined in fits that included the full correlation matrix, which was remade 
for each jackknife fit. The jackknife data sets on different ensembles were then combined 
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TABLE VII: Errors in the Kb, c parameters. The first column labels the heavy quark, the second 
gives the statistical and fitting error for the k parameter, the third gives the discretization error, 
and the fourth combines these in quadrature. 



K 


statistics + fitting 


discretization 


total 


K c 


1.2% 


0.3% 


1.2% 


Kb 


5.6% 


1.3% 


5.7% 



into a larger block-diagonal jackknife data set in order to perform the chiral fits. In this way, 
the fully correlated statistical errors were propagated through to the final result. 

With our high statistics (several hundred lattice gauge field configurations for each en- 
semble), we are able to resolve the full covariance matrix well enough that we do not need to 
apply a singular value decomposition cut on the eigenvalues of the covariance matrix. The 
double ratio fit is needed to establish h,A 1 (l) at the fiducial point (which was computed on 
the 0.0124/0.031 fine ensemble), while the quadruple ratios, lZ va i and lZ sca , are computed 
on the other ensembles in order to perform the chiral extrapolation and to remove taste 
breaking non-analytic terms. We find that the fit to the double ratio at the fiducial point on 
the 0.0124/0.031 ensemble is well described by a constant over a range of seven time slices. 
The excited state contamination in the quadruple ratios is even further suppressed, and we 
find that the correlated x 2 values allow for a constant fit region of six to ten time slices, 
depending upon the lattice spacing. We take the good correlated x 2 /d.o.f., ranging from 
0.15 to 1.00, in our constant plateau fits as evidence that the excited state contamination in 
these fits is negligible as compared to other errors. 

As an additional check of the jackknife fitting procedure, bootstrap fits were done to all 
of the double and quadruple ratios needed for this work. Close agreement was found for 
both central values and statistical errors. The statistical errors were typically the same size 
within 10%, and central values were well within la. The jackknife procedure had slightly 
larger errors than that of the bootstrap. 
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B. Heavy-quark mass dependence 



The value for depends on the heavy-quark masses, which are set by tuning the 

hopping parameters Kb and k c . The principal method starts by fitting the lattice pole energy 
to .E'(p) to the dispersion relation, 

2M 2 



E(p) = M 1 + ^- + b lP ' + b 2 J2W 4 + ■ ■■ , (40) 



in order to obtain the kinetic mass M 2 (as well as bi and 62, which are unimportant here). 



In the Fermilab method 



2n 



k is adjusted so that the kinetic mass agrees with 
experiment. Here we take the spin-average of kinetic masses of pseudoscalar and vector 
heavy-strange mesons and obtain our central values for Kb or k c , respectively, from the (spin- 
averaged B y s ' and Dy masses. Applying this procedure we find statistical and fitting errors 
of 5.6% for Kb and 1.2% for k c on the fine (a = 0.09 fm) ensembles. There is an additional 
error in k due to discretization effects. We determine this error by estimating the size of 
discretization effects for the Fermilab action (at a = 0.09 fm) as in Ref. [3]. This error is 
1.3% for Kb and 0.3% for k c . Adding in quadrature the statistical and fitting error together 
with the discretization error leads to a total relative uncertainty of 5.7% for Kb and 1.2% for 
k c . This error budget is summarized in Table IVHI Note that these errors are conservative 
and are likely to decrease substantially with more sophisticated fitting methods and the 
higher statistics data set currently being generated. 

We have computed /iAi(l) a t several different values of the bare charm and bottom quark 
masses, and these simulated points can be used to estimate the error in Ha x (1) from the above 
uncertainties in the tuning of the heavy-quark k values. Figure H] illustrates the dependence 
of ftvii(l) as a function of bottom and charm quark k values on one of the coarse (a = 0.12 
fm) ensembles. The points labelled Kb show fiA 1 {l) where we have fixed k c to the tuned 
charm value, but vary the bare Kb along the x-axis. The points labelled k c are similar, where 
the value of Kb is fixed at its tuned value, and the bare k c is varied. The above uncertainties 
in the k's, combined with the variation of /^(l) with k, lead to a systematic error of 0.7% 
in /iAi(l), labelled "kappa tuning" in Table IXl 
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FIG. 4: /i J 4 1 (l) for different values on the coarse fh' = 0.02 ensemble (full QCD point). The 
points labelled Kb show how /^(l) depends on when k c is fixed to its tuned value. For the 
points labelled k c the roles of and reversed. 



C. Perturbation theory 

The perturbative calculation of pA 1 is needed to match the heavy-quark lattice current, 
and the calculation has been carried out to one-loop order [0(a s )]. As discussed in Section 
IV, much of the renormalization cancels when forming the ratios of Z factors that define p 
[Eq. (I29p ], and the coefficients of the perturbation series are small, by construction. The one- 
loop correction is quite small, only 0.3-0.4% on the different lattice spacings. We take the 
entire one-loop correction of 0.3% on the fine lattices as an estimate of the error introduced 
by neglecting higher orders in the perturbative expansion. 
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D. Chiral extrapolation 



We estimate our systematic error due to the chiral extrapolation by comparing fits with 
and without the additional terms with coefficients q in Eq. fl33|) . i.e. analytic terms of higher 
order than NLO in rS^PT, since the two-loop NNLO logarithms are unknown. We also 
compare with continuum xPT, both NLO and (partial) NNLO. There are additional errors 
due to the uncertainties in the parameters that enter the NLO rS^PT formulas. By far the 
largest uncertainty of this kind is that due to the uncertainty in gDD*n- Finally, there is 
an error due to a mistuning of the parameter u on the coarse lattices. All of these errors 
are discussed below in more detail. In the discussion of chiral extrapolation errors, it is 
important to keep in mind that the chiral logarithms (either rS^PT or continuum) are tiny 
(~ 3 x 10~ 3 ) in the region where we have data. Non-analytic behavior is important only near 
the physical pion mass where the xPT should be a very good description in the continuum. 
The main feature of the chiral extrapolation is a cusp that appears close to the physical pion 
mass (in the valence sector), due to the Dn threshold and the fact that the D-D* splitting is 
very close to the physical pion mass. This cusp represents real physics, and must be included 
in any version of the chiral extrapolation used to estimate systematic errors. 

We extrapolate the light sea and light valence quark masses from the values used in the 
simulations, between m s /2 and m s /10, to the average physical light quark mass, around 
m s /27. We use staggered chiral perturbation theory and the prescription introduced in 
Section VI to remove the non-analytic taste-breaking discretization effects coming from the 
staggered light quark sector. Separate fits are performed for the two ratios introduced in 
Eqs. ( l36|) and ( l37l) . 1Z sea and 1Z va \. The chiral extrapolation is performed on these ratios, 
and the staggered discretization errors appearing in the NLO chiral logarithms are removed 
by taking a — > in the rS%PT expression. With the NNLO analytic terms given in Eq. ( 1331) 
the chiral extrapolation formulas for the ratios are 

ft va i = 1 + NLOi ogs + Cl m 2 Xp , (41) 
K se& = 1 + NL0 logs + c 2 {2m 2 Up +m 2 Sp ), (42) 

where NLOi ogs is a schematic notation representing the chiral logarithms coming from nu- 
merator and denominator. These terms are different for the two ratios, and can be obtained 
straightforwardly from the definitions of the ratios Eqs. (1361) and (I37|) . and the formula for 
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FIG. 5: 7£ V ai ratio versus valence pion mass squared on all ensembles for the three different lattice 
spacings. The curve is the continuum prediction through NLO in continuum xPT for this quantity. 
(See Appendix.) 



the non-analytic terms in Eq (lAip . The formula for 1Z W3 \ in the continuum is given explicitly 
in Eq. flA6[) . for the purposes of illustration. The NNLO term 03a 2 in Eq. (|33l) cancels in 
the ratios, and TZ sca , and TZ va \ each determine one of the remaining two NNLO coefficients. 
Note that the factor of tja in Eqs. ( I32I) and (!33l) cancels in the chiral formulas for the two 
ratios. The only free parameters in our chiral fits are c\ and Ci\ the rest are determined from 
phenomenology or from rS^PT fits to the pseudoscalar sector. 

The ratios in Eqs. (I3T)|) and (j3"7|) are completely predicted through NLO in the continuum 
once f n , Qdd*txi and the D-D* splitting are taken from experiment. The constants 

„2 

f n and gDD*ir appear in an overall multiplicative factor in front of the logarithmic 

term, as can be seen in Eq. flAlj) and Eq. (IA6I) . We take a fairly conservative range for the 
constant goD'-w determined from phenomenology, as discussed below, and the errors in this 
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quantity are accounted for in our final error budget. In the mass region where we have data, 
the NLO continuum chiral logarithms contribute to ^Ai(l) at the ~ 3 x 10~ 3 level or less. 
Figure [5] illustrates this, where the NLO continuum %PT prediction Eq. (1A6I) is plotted over 
our data points for 1Z V£L \. We find that the NLO continuum ;\TT describes the data quite 
well, giving a xVd.o.f. = 0.91 and a corresponding CL=0.51. This result is unchanged in 
the rS^PT fits; the effects of staggering are negligible in the region where we have data. 
We include the term proportional to C\ in Eq. (j4ip in our fits used to obtain the central 
value for this quantity, as explained in Section V. (Since including a linear term proportional 
to c\ increases the statistical error in h^ 1 , we take our central value and statistical error 
from this fit to be conservative.) This "partial NNLO" fit also has a good x 2 /d.o.f. = 1.05, 
with a corresponding CL=0.39. The constant linear term is small and consistent with zero 
[c\ = —0.006(15)]. Figure [6] shows the fit to 7Z vau \ versus m Xp for all three lattice spacings 
using the rS^PT formula, Eq. (JH|) . 

Although the data for 1Z va \ is consistent with a constant, the cusp appearing close to the 
physical pion mass is a prediction of NLO xPT and has a physical origin, namely the D-n 
threshold, as we have remarked. Thus, any fits used to estimate systematic errors, even 
those that are somewhat ad hoc, such as those including higher order polynomial terms, 
must include this cusp. Note that the cusp appears at the physical pion mass (in either 
SU(3) or SU(2) %PT), and is therefore in a region where xPT is expected to be a reliable 
expansion. The cusp is a property of the function F(m, /m) given in Eq. (1A2|) . and the 
position of the cusp as a function of m x is determined by the D-D* splitting A^ c ^ and the 
physical pion mass. We take these two quantities from experiment rather than from the 
lattice, since the experimental uncertainties are much smaller. 

We find that with or without the NNLO analytic terms, the %PT (continuum or rooted 
staggered) describes the lattice data with x 2 /d.o.f. close to 1 and correspondingly good 
confidence levels. We find a confidence level for the fit to 1Z sea of 0.76 for the fit that 
includes NNLO analytic terms. The strictly NLO expression for the lattice ratio TZ sca . has 
no free parameters, but it describes the data with a confidence level of 0.73. Similar fits to 
7?. va i are described above and yield reasonable confidence levels for both types of fits. Since 
the lattice data do not distinguish between these model fit functions, and the fit using only 
the NNLO analytic terms is not systematic in the chiral expansion, we assign the difference 
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FIG. 6: 72-vai ratio versus valence pion mass squared on all ensembles for the three different lattice 
spacings. The curve is the fit with 1 sigma error band to the ratio for all three lattice spacings using 
rSxPT, extrapolated to the continuum by taking a — > in the NLO staggered chiral logarithms. 

between the two determinations, which is 0.9%, as the systematic error of leaving out higher 
order terms when performing the chiral extrapolation. The final results for 7£ sea , TZ va \, and 
T^fid are given in Table lylTTl The errors are statistical only; note that the strictly NLO values 
have no free parameters, and therefore no statistical errors. The final value of still has 
statistical errors coming from the statistical errors in h^. The extrapolated results for TZm 
are consistent within the statistical errors of the NNLO fit. Again, we choose for our central 
value the result from the NNLO extrapolation with its larger errors to be conservative. 

The cyan (gray) band in Figure [6] is the continuum extrapolation with a — > in the rS^PT 
formula. For this quantity, the staggered lattice artifacts affecting the chiral logarithms in 
Ua x are negligible in the region where we have lattice data, which is due mainly to the small 
size of the chiral logarithms themselves. This is confirmed by the close agreement between 
the data points at each lattice spacing and the continuum curve. In fact, if we use continuum 
%PT to perform the chiral extrapolation, the result is unchanged. The primary difference 
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x extrapolated value 



TABLE VIII: Continuum extrapolated values of TZ SC& , 7£ va i, TZm, and evaluated at the 

physical quark masses. The first column labels the quantity. The second is the computed value 
including NNLO analytic terms in the chiral fit. The third is the quantity evaluated in purely NLO 
xPT, and has no free parameters (once gDD*%, fn an d A( c ) are taken from phenomenology) in the 
chiral fit. The final row shows ^ 1 (1), which includes a statistical error coming from h A d . The 
numbers are the same to the quoted precision using rS%PT or continuum xPT- 

w/ NNLO strictly NLO 



ftsea 1.0059(90) 0.9983 

7£ va i 0.9910(34) 0.9895 

Km 0.997(10) 0.9878 

h Al (l) 0.921(13) 0.9124(84) 



between the rS^PT expression and the continuum x?T expression is the reduction of the 
cusp near the physical pion mass in rSxPT, though our lattice data are not near enough to 
the physical pion mass to demonstrate this effect. 

Figure [7] shows the fit to lZ seau , extrapolated to the continuum and to the physical strange 
sea quark mass. Note that this ratio does not produce a cancellation of correlations between 
numerator and denominator and so has larger statistical errors than TZ va \. Here again the 
discretization effects due to staggered logarithms are negligibly small. Since the effects 
of including staggered discretization effects in the chiral logarithms are negligible in the 
region where we have numerical data, and since the only nontrivial feature in the chiral 
extrapolation is the cusp near the physical pion mass, which we describe by continuum x?T 
(our extrapolated curve has a — > in the rS^PT formula and thus reduces to the continuum 
form), we conclude that staggered taste- violating effects appearing in chiral logarithms are 
essentially removed in our ratio extrapolations. 

Figure [8] shows all of the full QCD points on the three lattice spacings. The curve is the 
quantity, 

hf^(m') « h%{mf,m hA ,mf ,a &d ) x [ft sca (m', m s , 0) x ft val (m', m', m s , 0)], (43) 
which again becomes an exact relation for the physical form factor when a — > 0. The curve 
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FIG. 7: 7£ sca ratio versus mfj p for all ensembles and lattice spacings. The curve is the fit to all 
of the lattice data, extrapolated to the continuum. The curve is also extrapolated to the physical 
strange sea quark mass. 



is thus the product of the two continuum extrapolated ratio fits shown in Figures [6] and [7J 
times the fiducial point, which we take to be am fid = 0.0124 at the fine lattice spacing (the 
solid square in Figure Because this is a full QCD curve, the valence mass m x equals 
the light sea mass m! . The other full QCD points are shown as open symbols in Figure [8] 
for comparison, though the fits were performed on the ratios and normalized by the fiducial 
point at am fid = 0.0124. Note that the curve is already extrapolated in the strange sea 
quark mass, and so does not perfectly overlap with the om M = 0.0124 point. As discussed 
above, when this quantity is evaluated at fh' = m it yields the value of Ha x at physical quark 
masses. The cross is the extrapolated value, where the solid line is the statistical error, and 
the dashed line is the total systematic error added to the statistical error in quadrature. 

The low energy constant gDD*n enters the chiral extrapolation formula and determines 
the size of the cusp near the physical pion mass. Our data do not constrain this constant, 
so we take a wide range for goD*n that encompasses the range of values coming from phe- 
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m n 2 (GeV 2 ) 

FIG. 8: The full QCD points versus m% on the three lattice spacings are shown in comparison 
to the continuum curve. The curve is the product of the two continuum extrapolated ratio fits 
shown in Figs. (|6|) and (J7J), times the fiducial point, which we have chosen to be the fh! = 0.0124 
fine lattice point (the filled square). The curve is already extrapolated to the physical strange sea 
quark mass, and so does not perfectly overlap with the lattice data point at the fiducial value. The 
cross is the extrapolated value, where the solid line is the statistical error, and the dashed line is 
the total systematic error added to the statistical error in quadrature. 



nomenology and lattice calculations: fits to a wide range of experimental data prior to the 



n 97+O.O6 
u - z ' -0.03 



measurement of the D* width by Stewart (goD* 
art analysis including the D* width (goD*-* = 0.51; no error quoted_ 
(9 



73|), an update of the Stew- 



N f =0 



c 0.38 [75j), quenched lattice QCD (q p p ,^ - 

N-—^ ' b 



74J), quark models 



0.67±0.08±He M, two flavor lattice 



QCD in the static lim it ((/.,.;.:,. = 0.516 ± 0.051 (77|), and the measurement of the D* width 
(gDD*ir = 0.59 ± 0.07 [78j). There are as of yet no 2+1 flavor lattice calculations of gDD*n- 



For this work we take gr>D*ix — 0.51 ± 0.2, leading to a parametric uncertainty of 0.9% in 
/iAi(l) that is included as a systematic error. 
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The additional low energy constants that enter the chiral formulas are the tree-level 
continuum coefficients /io and /, and the taste- violating parameters that vanish in the con- 
tinuum. These are the taste splittings, a 2 An with S = P, A, T, V, I, and the taste- violating 
hairpin-coefficients, a 2 8' A and a 2 5' v . We set / to the experimental value of the pion decay 
constant, f n = 0.1307 GeV, in the coefficient of the NLO logarithms. The pion masses used 
as inputs in the rS^PT formulas are computed from the bare quark masses and converted 
into physical units using 

m ly = ( r 1 / r l hyS ) 2 three (m x +m y ), (44) 

where /i tree is obtained from fits to the light pseudo-scalar mass squared to the tree-level form 
(in n units), rlfi tree (m x +m y ). This accounts for higher-order chiral corrections and is more 
accurate than using \i obtained in the chiral limit, giving a better approximation to the pion 
mass squared at a given light quark mass. Since the parameters in our lattice simulations 
at different lattice spacings are expressed in t\ units, we require the physical value of t\ to 
convert to physical units and take the physical pion mass and from experiment. Thus, 
the ~ 2.5% uncertainty in r^ hys gives a parametric error in the chiral extrapolation. Because 
the chiral extrapolation is so mild, however, this error turns out to be negligible compared 
to other systematic errors. Since we are taking the pion mass from experiment there is 
a negligible error due to the light quark mass uncertainty in the chiral extrapolation. The 
strange sea quark mass enters the chiral extrapolation formulas, but the dependence is weak, 
and the error in the bare strange quark mass leads to a negligible parametric error in Jia 1 - 



The taste-splittings A= have been determined in Ref. 16], and their approximately 10% 
uncertainty also leads to a negligible error in /^(l). The taste-violating hairpin coefficients 
have much larger fractional uncertainties, but these too lead to a negligible uncertainty in 
tiA 1 (l)- Even setting the rS^PT parameters to zero does not change our result for h Al (l) 
significantly. As mentioned above, our result does not change if we use the continuum x?T 
formula in our chiral fits. 

In the calculation of the form factor, the tadpole improved coefficient Csw = l/ M o ^ s 
obtained with u$ from the Landau link on the coarse lattices, but from the plaquette for 
Mo on the fine and coarser lattices. Though unintentional, there is nothing wrong with this, 
since it is not known a priori which provides the best estimate of the tadpole improvement 
factor. However, the Uq term for the spectator light (staggered) quark, which appears in the 
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tadpole improvement of the Asqtad action, was taken from the Landau link on the coarse 
lattices, even though the sea quark sector used uq from the plaquette. On the fine and 
coarser lattices, uq was taken to be the same in the light valence and sea quark sectors. The 
estimates of u from plaquette versus Landau link differ only by 4% on the coarse lattices. 

Although the effect of this mistuning is expected to be small (correcting no would lead to 
a slightly different valence propagator and different tuned k values, thus leading to a small 
modification of the staggered chiral parameters in the valence sector for the coarse lattices 
used as inputs to the chiral fit), it is possible to study how much difference it makes using 
the Jia-l lattice data. Including all three lattice spacings and using our preferred chiral fit, 
we find h^^l) = 0.921(13) where the error here is statistical only. If we neglect the coarse 
data points, we find = 0.920(17), almost unchanged except for a somewhat larger 

statistical error. We can also examine the ratios TZ va \ and lZ se& . In our preferred fit to all 
the lattice data these are 0.9910(34) and 1.0059(90) respectively, where the errors are again 
only statistical. If we drop the coarse lattice data, these become 0.9960(56) and 0.999(13) 
respectively. Since the ratio TZ sea has very little valence quark mass dependence, we can 
combine 7£ sea from the fit to all of the lattice data with TZ VSl i from the fit neglecting the 
coarse lattice data. This is useful, because TZ sea , has the larger statistical error, so we would 
like to use the full lattice data set to determine this ratio, thus isolating the mistuning in 
the valence sector on the coarse lattices. When this is done we find that the central value 
of the final h^^l) is shifted upward by 0.4%, well within statistical errors and smaller than 
our other systematic errors. We assign a systematic error of 0.4% due to the Uq mistuning. 



E. Finite volume effects 



The finite volume corrections to the integrals which appear in heavy-light xPT formulas 



including those for B — > D* were given by Arndt and Lin [79J. There are no new integrals 
appearing in the staggered case, and it is straightforward to use the results of Arndt and 
Lin in the rS^PT for /^(l), as shown in Ref. [22(]. We find that although the finite volume 
corrections in /lAi(l) would be large near the cusp at the physical pion mass on the current 
MILC ensembles (ranging in size from 2.5-3.5 fm), for the less chiral data points at which we 
have actually simulated, the finite volume effects are negligible. For all data points in our 
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TABLE IX: hA l (l) at physical quark masses at different lattice spacings, where taste- violating 
effects have been removed, or shown to be negligible. Discretization effects due to analytic terms 
associated with the light quark sector and heavy-quark discretization effects remain in the lattice 
data. 



a (fm) 


h Al (l) 


0.15 


0.914(11) 


0.12 


0.907(14) 


0.09 


0.921(13) 



simulations the finite volume corrections are less than 1 part in 10 4 . We therefore assign no 
error due to finite volume effects. 

F. Discretization errors 



As shown in Ref. 
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23, 



ad, 



49], the matching of lattice gauge theory to QCD is accom- 
plished by normalizing the first few terms in the heavy-quark expansion. This is done by 
tuning the kinetic masses of the D s and B s mesons computed using the SW action (for the 
heavy quarks) to the experimental meson masses. Tree-level tadpole- improved perturbation 
theory is used to tune the coupling csw and the rotation coefficient d\ for the bottom and 
charm quarks. Once this matching is done, the discretization errors in h^iV) are of order 



a s (A/2mQ) 2 and (A/2mQ) 3 28j], where the powers of two are combinatoric factors. The 
leading matching uncertainty is of the order a s (A/2m c ) 2 . We estimate the size of this error 
setting a s = 0.3, A = 500 MeV, and m c = 1.2 GeV, which gives a s (K/2m c f = 0.013. 

Since we have numerical data at three lattice spacings we are able to study how well 
the power counting estimate accounts for observed discretization effects. Making use of 
Eq. (l4"3l) . but varying the fiducial lattice spacing from our lightest to coarsest lattices, we are 
able to obtain (1) at physical quark masses, with discretization effects associated with the 
staggered chiral logarithms removed in the ratios appearing in Eq. (j4"3"|) . The discretization 
effects that remain are: taste- violations in , taste violations at higher order than NLO 
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FIG. 9: /ia^I) at physical quark masses versus a? (fm 2 ) where taste-violating effects have been 
removed, or shown to be negligible. Discretization effects due to analytic terms associated with the 
light quark sector and heavy-quark discretization effects remain in the lattice data. 



in the ratios, the effect of the analytic term coming from light quark discretization effects 
(proportional to a s a 2 ), and the heavy-quark discretization effects. The taste- violations in 
and the taste-violations in the ratios appearing at higher order than NLO have been 
shown to be negligible. We now consider the remaining discretization errors coming from the 
light quark analytic term and the heavy-quark discretization effects. Table IIXI presents the 
results for /^(l) as obtained from Eq. (]43j) and Figure [9] shows them plotted as a function 
of lattice spacing squared. Although the Fermilab action and currents possess a smooth 
continuum limit, the MILC ensembles are not yet at small enough a to obtain simply 0(a) 
or 0(a 2 ) behavior. The spread of the lattice data points gives some indication of the size 
of the remaining discretization effects, however, and we find that the fine (0.09 fm) lattice 
data point and the coarse (0.12 fm) lattice data point differ by 1.5%. This is similar to our 
power counting estimate, and we assign the larger of the two, 1.5%, as the systematic error 
due to residual discretization effects. 
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G. Summary 



Our final result, given the error budget in Table IXl is 

h Al = 0.921(13)(8)(8)(14)(6)(3)(4), (45) 

where the errors are statistical, parametric uncertainty in gDD*n, chiral extrapolation errors, 
discretization errors, parametric uncertainty in heavy-quark masses (kappa tuning), pertur- 
bative matching, and the uq (mis)tuning on the coarse lattices. Adding all systematic errors 
in quadrature, we obtain 

h Al (l) = 0.921(13) (20). (46) 



This final result differs slightly from that presented at Lattice 2007 



80j , where a preliminary 



1ia 1 {1) = 0.924(12)(19) was quoted. There are three main changes in the analysis from 
the preliminary result: our earlier result used a value of a s in the perturbative matching 
evaluated at the scale 2/a, while the present result uses the HLM 64j prescription to fix the 
scale. This causes a change of 0.1%, well within the estimated systematic error due to the 
perturbative matching. In the previous result, the fine lattice data was blocked by 4 in the 
jackknife procedure; we now block by 8 to fully account for autocorrelation errors. This does 
not change the central value, but increases the statistical error slightly. Finally, we have 
chosen a value for gnD*n = 0.51 ± 0.2 instead of gDD*n = 0.45 ± 0.15 to be more consistent 
with the range of values quoted in the literature. This causes a decrease in h Al (l) of 0.2%. 

VIII. CONCLUSIONS 

We have introduced a new method to calculate the zero-recoil form factor for the B — > 
D*lv decay. We include 2+1 flavors of sea quarks in the generation of the gauge ensembles, 
so the calculation is completely unquenched. We have introduced a new double ratio, which 
gives the form factor directly, and leads to a large savings in the computational cost. The 
simulation is performed in a regime where we expect rooted staggered chiral perturbation 
theory to apply; we therefore use the rS^PT result for the B — > D* form factor |22j] to perform 
the chiral extrapolation and to remove taste-breaking effects. To aid the chiral and continuum 
extrapolations, we introduced a set of ratios that has allowed us to largely disentangle light 
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TABLE X: Final error budget for hA 1 (l) where each error is discussed in the text. Systematic 
errors are added in quadrature and combined in quadrature with the statistical error to obtain the 
total error. 



Uncertainty 


^i(l) 


Statistics 


1.4% 


9DD*tt 


0.9% 


NLO vs NNLO xPT fits 


0.9% 


Discretization errors 


1.5% 


Kappa tuning 


0.7% 


Perturbation theory 


0.3% 


uq tuning 


0.4% 


Total 


2.6% 



and heavy-quark discretization effects. Our new result, ^(l) = Ha, (1) = 0.921(13) (20) is 
consistent with the previous quenched result, .F(l) = 0.913±g;g§| [13] , but our errors are both 
smaller and under better theoretical control. This result allows us to extract \V c b\ from the 
experimental measurement of the B — > D*iv form factor, which determines ^(1)1 V^l- After 
applying a 0.7% electromagnetic correction to our value for ^(l) 81[, and taking the most 
recent PDG average for \V ch \T{\) = (35.9 ± 0.8) x 10" 3 [SsJ, we find 



\V cb \ = (38.7 ±0.9 cxp ± 1.0 



■"thco , 



x 10" 



(47) 



This differs by about 2a from the inclusive determination \V c b\ = (41.6 ± 0.6) x 10 3 [82| 
Our new value supersedes the previous Fermilab quenched number 13j, as it should other 



quenched numbers such as that in Ref. 



Our largest error in .F(l) is the systematic error due to heavy-quark discretization effects, 
which we have estimated using HQET power counting and inspection of the numerical data 
at three lattice spacings. This error can be reduced by going to finer lattice spacings, or 
by using an improved Fermilab action 70|]. When using this improved action, it would be 



calculates the B — > D*tv form factor in the quenched approximation at zero and non-zero recoil 



Ref. 

momentum and uses a step-scaling method 8J] to control the heavy-quark discretization errors. 



41 



necessary to improve the currents to the same order. We have introduced a method for 
separating the heavy and light-quark discretization errors, where the physical Ha x can be 
factorized into two factors, h A d x lZ&d, such that the heavy quark discretization errors are 
largely isolated in h A d . Combining our value of TZ^d = 0.997(10) (13) (where the first error 
is statistical, and the second is due to systematics that do not cancel in the ratio) with a 
determination of h A d at finer lattice spacings and/or with an improved action would be a 
cost-effective way of reducing the heavy-quark discretization errors. The next largest error 
in our calculation of is statistical, and this error drives many of the systematic errors. 
This is mostly a matter of computing. It would also be desirable to perform the matching of 
the heavy-quark current to higher order in perturbation theory, or by using non-perturbative 
matching. With these improvements, it would be possible to bring the error in !F(1) to or 
below 1%, allowing a very precise determination of \V c t,\ from exclusive semi-leptonic decays. 

Acknowledgments 

We thank Jon Bailey for a careful reading of the manuscript. Computations for this work 
were carried out in part on facilities of the USQCD Collaboration, which are funded by the 
Office of Science of the U.S. Department of Energy; and on facilities of the NSF Teragrid un- 
der allocation TG-MCA93S002. This work was supported in part by the United States 
Department of Energy under Grant Nos. DE-FC02-06ER41446 (CD., L.L.), DE-FG02- 
91ER40661(S.C), DE-FG02-91ER40677 (A.X.K.), DE-FG02-91ER40628 (C.B., J.L.), DE- 
FG02-04ER41298 (D.T.) and by the National Science Foundation under Grant Nos. PHY- 
0555243, PHY-0757333, PHY-0703296 (CD., L.L.), PHY-0555235 (J.L.), PHY-0456556 
(R.S.). R.T.E. and E.G. thank Fermilab and URA for their hospitality. Fermilab is op- 
erated by Fermi Research Alliance, LLC, under Contract No. DE-AC02-07CH11359 with 
the United States Department of Energy. 



42 



APPENDIX A: CHIRAL PERTURBATION THEORY 



Eq. (34) of Ref. |22| gives the expression needed for h Al (l) in partially-quenched 
with degenerate up and down quark masses (the 2+1 case) in the rooted staggered theory: 

9 2 DD*n I 1 



(B X )PQ,2+1 



1)/ VA = 1+X A (A X ) + 



E ^ 



48vr 2 P 16 

j=xu,xu,xs 
Z=I ,P,4V,4A,6T 



J6{M} 5) } 



+ a 2 5' v 



{M) 



a. 



2€{M<, 7) } 



£ bK,({m«};{M)F, 

(Al) 



where 



F K, = -1 <j ^ In + - Az 3 + 2tt 



*?-l(z? + 2)fln 



1 - 2%(^- 



$-1) 



(A^) 2 ln(^)+0[(A( c )) 3 ], 

X • 



(A2) 



with F(rrij, Zj) = F(rrij, —Zj), and Zj = A^°' /mj, where is the D-D* mass splitting. The 

19| . and for completeness we quote them 



residues i?[ n '^ and D^ 1 -'^ are defined in Refs. 
here: 



R [ ?' k] ({M},{n}) 



, 2 y 



dm 2 j 



D [ -i k] ({M},{ri) 



M n,fcl (W,M). 



(A3) 



These residues are a function of two sets of masses, the numerator masses, {M} = 
{m 1; m 2 , m n } and the denominator masses, {/i} = {/ii, /12, fik}- in our 2+1 flavor 



case, we have 



{M| } } = {m„ mx }, 
{M^ 7) } ee {m v ,m v ,,m x }, 
{/j,} ee {m[/,m 5 }. 



(A4) 
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The masses m m , m Vv , m v i are given by 



18| 
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The ratio in the continuum through NLO in x?T is 
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where 



(A5) 



(A6) 



(AT) 



and where my is a valence pion made of two quarks set to the fiducial valence quark mass, 
and the subscript x' refers to a valence quark at the fiducial mass. This ratio is one by 
construction when the valence quark mass equals the fiducial valence quark mass. 
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